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We construct in the K matrix formalism concrete examples of symmetry enriched topological 
phases, namely intrinsically topological phases with global symmetries. We focus on the Abelian and 
non-chiral topological phases and demonstrate by our examples how the interplay between the global 
symmetry and the fusion algebra of the anyons of a topologically ordered system determines the 
existence of gapless edge modes protected by the symmetry and that a (quasi)-group structure can 
be defined among these phases. Our examples include phases that display charge fractionalization 
and more exotic non-local anyon exchange under global symmetry that correspond to general group 
extensions of the global symmetry group. 
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I. INTRODUCTION 

The understanding of phases of matter has come a long 
way beyond the Landau paradigm-'^. Different phases of 
matter cannot be simply classified by Landau's symme- 
try breaking and a corresponding order parameter. In 
the case of phases involving short-range entanglement 
(SRE), it is now realized that for a given preserved 
global symmetry G^, they could be subdivided into many 
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different phases that cannot be connected by any lo- 
cal, unitary transformation without breaking the sym- 
metry. These phases, called Symmetry Protected Topo- 
logical (SPT) phases, turn out to be classified by group 
cohomology H^{Gs,U{l)) in two spacetime dimensions 
and are believed to be classified by H'^{Gs,U{l)) in d- 
dimensional spacetime^. An independent study based 
on K matrix construction that is particularly powerful 
in studying Abclian symmetry groups have confirmed 
many of the group cohomology classification, and more- 
over shed light on the edge excitations and transport 
properties of these phases^'^. Things are more interest- 
ing when phases possess long range entanglement (LRE). 
Even without symmetry, they already show a very rich 
bulk structure, and so far only a partial classification of 
them is known. It is expected that when symmetry is 
incorporated, where such phases arc often dubbed Sym- 
metry Enriched Topological (SET) phases, they would be 
further subdivided into different phases, and the allowed 
action of the symmetry group can be very exotic as it is 
already anticipated in earlier work on projective symme- 
try group where charge fractionalization is one common 
feature there^. 

More recently, there is renewed interest in systemati- 
cally constructing and classifying these SET phases, no- 
tably in Ref^"^ . Here, we would like to extend the meth- 
ods in Ref^'^ to constructing SET phases in 2-1-1 dimen- 
sions. Our extension confirms many of the results in 
Rcf^^*, particularly regarding the conditions of charge 
fractionalization and generalization to non-local symme- 
try transformations. Moreover, the K matrix analy- 
sis allows us to study the edge excitations in the pres- 
ence of boundaries. A (quasi)-group structure among 
phases with the same global symmetry and fusion alge- 
bra emerges as we consider stacking them together, which 
does not appear to be directly related to group cohomol- 
ogy, although such a relation was found in the case of 
SPT phases'^. We also generalize constructions in Ref"^ 
to include some non-Abelian symmetry groups. 

The K matrix construction is most powerful in deal- 
ing with phases whose anyons are governed by Abelian 
statistics. However, our study has inspired us of a more 
general way to construct and perhaps ultimately to clas- 
sify symmetric phases with non-Abelian anyons. We will 
comment on the general idea towards the end of the pa- 
per. 

Since our construction amalgamates and generalizes 
several ideas, we would like to begin our discussion with 
a general overview that puts together the various build- 
ing blocks necessary for the current paper and clarify a 
few core concepts. 

In section HI we will introduce our K matrix approach, 
based on the ideas developed in Ref^'"' how K matrix 
can be taken as the starting point for constructing SPT 
phases. Then we apply this approach in section IV to 
construct LRE phases with symmetries. We will dwell 
particularly on the symmetry enriched Z2 gauge theory 
and the double semion model, studying their edge excita- 



tions and a quasi-group structure that emerges between 
the phases. This is then generalized in section V to Zjv/ 
symmetry in phases with Z^r x Zyv and related fusion 
algebras. 

To explore more exotic group actions of the symmetry 
group beyond charge fractionalization, we study some ex- 
amples in VI that involve anyon exchange based again on 
phases with Z2 x Z2 fusion algebra. More exotic and elab- 
orate examples based on phases with fusion Z4 x Z4 is 
discussed in section VII. We collect these ideas, and sum- 
marise the unifying principles behind these examples in 
section VIII, where we construct also new phases accom- 
modating discrete non-Abelian group actions, explicitly 
the Dihedral groups. 

We compare our results with existing results in the 
literature in section IX, and then conclude our discussion 
in section X with open questions. 

Appendices A and B collect some technicalities. Ap- 
pendix C, however, provides the K matrix version of the 
"duality" relation between a SPT phase and a topological 
phase, where the latter descends from gauging the global 
symmetry in the former. This relation was first proposed 
and realized in the string-net formalism in Ref^*' and then 
further discussed in Ref^'*'^^. 



II. SYMMETRY ENRICHED PHASES IN 2 + 1 
DIMENSION: AN OVERVIEW 

The main focus of the current paper is to construct 
examples where topological phases - namely phases that 
possess LRE - are endowed with global symmetries. The 
theme has received much attention recently, for example 
in Ref^^*. Several principles underlie these discussions 
and constructions, and we would like to summarize them, 
along with a conceptual account for our approach, before 
moving on to our explicit constructions that concretely 
realize these principles. 

Any discussion of symmetry can hardly avoid the in- 
troduction of groups. Since a number of groups would 
be introduced in our discussion, we would like to cata- 
log them here and explain briefly the role each plays, for 
clarity and later convenience. 

One important feature of LRE-phases is the emer- 
gence of non-local deconfined quasi-particlcs. In 2-1-1 di- 
mensions for instance, quasi-particles (anyons) displaying 
Abelian or non-Abelian anyonic statistics furnish such 
examples. While physical observables are characterized 
by local, bosonic excitations, anyons are non-local and 
cannot be physically excited in complete isolation. Nev- 
ertheless, the phases often exhibit "deconfined" limits, in 
which it is possible to keep various anyons so far apart 
that a lot of the operators can be considered as acting lo- 
cally on individual anyons present. This is analogous to 
the familiar situation in gauge theories, in which physical 
excitations are necessarily gauge invariant, even though 
it is often useful to think of them as composites being 
made up of gauge charged particles, particularly in a 
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"deconfined" limit when the charged particles can be, 
to some extent, isolated. In fact, a lot of these LRE 
phases can be conveniently described by gauge theories, 
such as the familiar case of Z2 spin-liquid, where a Z2 
gauge symmetry effectively emerges in the "deconfined" 
limit. Moreover, while it is unclear whether a complete 
classification of these LRE phases exist, and very likely, 
any such complete classification would invoke the Math- 
ematics of tensor categories^^'^^, the framework of gauge 
theories alone already encompasses a large class of LRE- 
phase^'^^'"'^^^^^ including many of the well-known paradig- 
matic examples such as the Z2 spin liquid. Most of the 
examples discussed in this paper are within the gauge 
theory framework, and thus wc will frequently refer to 
"gauge groups" Q in this sense. These gauge theories will 
be taken as the starting point on which we impose global 
symmetries. This starting point enables us to charac- 
terize or label an anyon- each topological sector- by its 
gauge charge and flux. A flux is labeled by a conjugacy 
class of Q, while the associated charge takes value in the 
irreducible (projective) representations of the centralizer 
subgroup of the flux in Q. Each anyonic excitation for 
any given gauge group would fall into one of the three 
categories: pure charge, pure flux, and dyon, which 
has both a flux label and the associated charge label. In 
the rest of the paper, we refer to different anyons using 
these terminologies where appropriate. Actually, as one 
will see, when a global symmetry is incorporated, an- 
other group that behaves essentially like a gauge group 
may appear, as wc will explore shortly below. This group 
however is generally different from Q. 

As discussed in the previous paragraph, 2 + 1 dimen- 
sions LRE phases generally bear anyonic, low-energy ex- 
citations, Abelian and/or non-Abelian. The interactions 
of these anyons are described by a set of fusion rules, in 
the sense that when viewed from sufficiently far away, 
various anyons relatively close together can be treated as 
a single lump. The lump behaves essentially as another 
anyon, now with a different topological charge and/or flux 
that descend from those of the constituent anyons in the 
lump. These fusion rules generally form an algebra ^, a 
fusion algebra, which in the case of Abelian anyons, is 
in fact an Abelian group. As discussed in Ref^ and later 
generalized in Ref® , ^ plays a central role in determining 
the possible ways a global symmetry could act. In this 
paper, like in most other discussions of symmetries, the 
global symmetries form a group Gg ■ Since symmetry acts 
reasonably locally in many cases, they can be understood 
as acting on individual anyons roughly independently. As 
emphasized above, however, anyons are not physical ex- 
citations and arc thus not directly a physical observable; 
therefore, it is conceivable that the physical states must 
transform linearly under the global symmetry - particu- 
larly that means they must transform trivially under the 
identity operator of the symmetry- such a restriction can 
be lifted on individual anyons. The simplest possibility 
is that the anyons live in projective representation spaces 
of Gs, in which case the anyons are considered to have 



undergone charge fractionalization. 

There are more exotic likelihoods, as demonstrated in 
Ref^'^ and also in some of our examples constructed in 
this paper in sections VI and VII, where exchange of 
anyons are involved and such symmetry transformations 
are not strictly local as opposed to the case of fraction- 
alization. Nonetheless, in all these cases, the fusion alge- 
bra/group 5" constrains admissible ways the anyons can 
transform, by demanding that the aggregate transfor- 
mations on any group of anyons that fuse to a physical 
bosonic excitation must be reduced to those correspond- 
ing to a linear representation, such that the identity op- 
erator acts trivially. In other words, we are effectively 
"modding out" transformations on anyons when the ag- 
gregate transformation of the group of anyons that fuse to 
a boson is trivial. These transformations that are modded 
out constitute a linear representation of a subgroup Ng 
of 5. In this sense therefore, Ng also behaves very much 
like some kind of gauge group, although it should not be 
confused with Q introduced earlier. They are generally 
different. 

To concisely describe and thus classify these non-trivial 
transformations, we can introduce yet another group G. 
In Ref^, G is the central extension of Gg by Ng. In that 
case, elements of Ng necessarily commute with those of 
Gs ■ This has been generalized to other group extensions, 
where Gs is the quotient subgroup G/Ng, where Ng is 
the normal subgroup of G. Anyons transform as linear 
representations of G, and these group extensions pro- 
vide the platform of classifying projective, and actually 
more general non-linear representations of Gs in which 
the anyons may fall into. In this fashion, the group ac- 
tions even of an Abelian Gs do not necessarily commute, 
examples of which have been seen in Ref^ and will be 
shown in this paper. And more generally, the group G 
is itself non-Abelian, and we obtain, to our knowledge, 
the first of such an example implementing non-Abelian 
group action in the K matrix construction, as discussed 
in section VIII. 

It should be noted that such classification of physi- 
cally admissible non-trivial actions of global symmetries 
Gs on any non-local excitations have appeared elsewhere. 
Most notably, in fermionic symmetry protected topolog- 
ical (SPT) phases , which involve only short-range en- 
tanglement, fermions nevertheless can transform projec- 
tivcly under Gg as long as any pair of them transform 
linearly. Framing it in the language developed above, the 
fusion group can be taken as ^5 = Z2 and the projective 
representations can be understood as group extension of 
Gs by Z2. In bosonic SPT phases, since the underlying 
excitations are already physical bosons, there is no no- 
tion of charge fractionalization, and in our language, the 
fusion group can be thought as = Zi, the trivial group. 

Before we close our discussions on Abelian phases, let 
us comment that the classification of symmetric LRE 
phases via the idea of group extensions does not a pri- 
ori inform us whether a given phase possesses non-trivial 
edge excitations in the presence of a boundary. Here, 



4 



non-trivial edge excitations refer to the edge modes of 
the anyons that cannot be gapped out without break- 
ing the symmetry and thus remain gapless as protected 
by the symmetry. The virtue of an expHcit construc- 
tion using K matrices is that one is able to explore the 
fate of the edge states as much as classifying them. De- 
spite transforming in highly non-linear representations 
under G^, there is no guarantee that the edge behaves 
also non-trivially. We found examples in which even very 
exotic transformation rules, implying charge fractional- 
ization and more, can lead to a gapped edge that respects 
the global symmetry. Among those phases that do pos- 
sess non-trivial edges, which feature gapless excitations 
or spontaneousely broken global symmetry Gg , it appears 
that there is a notion of a (quasi)-group structure be- 
tween them, when one considers stacking them together 
as in Rcf"^. This is a quasi-group also because the iden- 
tity is not a single clement but contains those phases 
that have fully gapped edge state without breaking Gg- 
This is discussed in section IV C. It is yet not completely 
clear whether such a group structure can always appear 
for any gauge theories, or that they are related to group 
cohomology, as in the case of SPT phases'^. 

In acknowledging the central and similar role ^ plays 
in constraining admissible non-trivial Gg representations 
in phases both short range and long range entangled, and 
also the quasi-group structure that ties together several 
phases that share the same fusion algebra we deem it 
convenient to refer generally to these phases as Symmetry 
Enriched Phases (SEP) and label classes of them with the 
same symmetry group Gg and fusion algebra ^ that are 
related by the group structure of SEP{^,Gg). This is 
in fact a unified notion of phases with symmetry that 
also encompass SPT phases: Fermionic SPT phases are 
classified by SEP{Z2,Gg), and bosonic SPT phases by 
SEP{Zi,Gg). 

Finally, let us comment on the situation of non-Abelian 
anyons. In the above discussion, we have very much re- 
stricted our attention almost entirely to Abelian anyons, 
whose fusion 5 is an Abelian group. This is also the ma- 
jor focus of our paper, where we make heavy use of the 
K matrix construction, which is appropriate for Abelian 
anyons. But, the discussion here, and also the discus- 
sion of group extensions discussed in Ref® have pointed 
to a general way to constructing LRE phases with sym- 
metries, if not completely classifying them. The idea is 
very much like the case of Abelian anyons, where different 
phases can be thought of as different embedding of the 
fusion group ^ inside a larger group G. Quite generally, 
particularly in the framework provided in Ref^'^^'^^"^^ 
describing large classes of LRE phases where the fusion 
^ forms a representation ring of a quantum group, which 
is an algebra it, an LRE phase possessing symmetries can 
be thought of as embedding il within a larger quantum 
group HQ. Analogous to the case of Abelian anyons, the 
quotient algebra is then taken as the global symmetry. 
This framework provides a natural way in which anyons, 
which fall into irreducible representations of the larger al- 



gebra, can be decomposed as a direct sum of irreducible 
representations of il, which in turn dictates how anyons 
transform under the global symmetry given by the quo- 
tient algebra. The embedding of a smaller invariant sub- 
group in a larger one employs the same Mathematics as 
in symmetry breaking, in which a large (gauge) group is 
broken to its invariant subgroup. For non-Abelian topo- 
logical phases, the relevant Mathematics would be that 
employed in Hopf symmetry breaking, which has been 
discussed in Ref^^"^° in the context of anyon condensa- 
tion. Many ideas can be directly applied here. We shall 
report a more detailed discussion elsewhere^ ^. an 



III. SYMMETRY ENRICHED PHASES IN 2 + 1 
DIMENSION: THE APPROACH 

In this section, we shall elaborate on our approach 
for studying LRE Abelian phases with symmetry. We 
take the formalism known as K matrix plus Higgs terms. 
This formalism was used in Ref'^ for studying bosonic and 
fermionic SPT phases in 2 -I- 1 dimension. We shall first 
briefly review the relevant pieces of this formalism then 
extend it to the case of LRE phases with symmetry. 



A. The K matrix + Higgs term formulation 

It is believed that 2+1 dimensional Abelian topological 
phases, including SRE phases and LRE Abelian, can be 
described in a unified fashion as effective Chern-Simons 
(CS) theories in the ii'-matrix formulation due to Wen et 
^^22-25 ^ -^^jjose generic Lagrangian density reads 

^cs = -^alKjjd^aje^'^P - aj^j^ • • • , (1) 

where Einstein's summation rule is assumed. The inter- 
nal indices I, J = 1,2,. ..,A^ label a set of N internal 
U{1) gauge fields a^, where the greek letters are space- 
time indices. The K matrix satisfies Kjj = Kjj S Z. 
A generic quasiparticle, however, is a fusion of the fun- 
damental ones and may be characterized by an integer 
vector 1 = {li, I2, ■ ■ ■ , In)'^ , carrying // unit of charge. 
The self statistics of a quasiparticle and the mutual statis- 
tics of two different quasiparticles a and b are respectively 
given by 

A physical quasiparticle is a boson, characterized by 
a vector 1^ that satisfies ^b/tt = and 6*50 /tt = 
(mod 27r) with arbitrary quasiparticle la. The ground 
state degeneracy (GSD) of the system placed on a torus 
is given by 

GSD = |A"|, (3) 
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where as an abuse of notation, K is the determinant of 
the K matrix. The Lagrangian in Eq. (1) describes SRE 
phases if = 1 and LRE phases if |X| > 1. In this 
paper, we concentrate on the latter case. Moreover, if 
K has the same number of positive and negative eigen- 
values, which also implies dim A' e 2Z+, it describes a 
non-chiral topological order. 

In the absence of any symmetry, one can condense the 
bosons by adding to Ccs potential terms: 

C'cs^Ccs+ E {c,\{b\' +h.c\ (4) 

l^bosons I 

where each C\ is constant, and 6/ is the annihilation op- 
erator of the fundamental excitation of type, with 
b\ = bj^. Each such term C\Y[ib\' -I- h.c. is often called 
a Higgs term. Note that this addition does not affect any 
topological properties of the system described by Cos- 
The K matrix theory makes it handy to study the edge 
states if the system has a boundary. The effective action 
of the edge theory is given by 

Se^S% + Si,, (5) 

where 

j AtAxY,{Kijdt4>id^4>J~Vijd^4>id^4>j) (6) 

corresponding to Ccs is the effective description of the 
gapless edge excitations, with 0/ the edge field associated 
with and Vjj a constant, positive definite matrix that 
determines the velocity of the edge excitations, and 

^E= C{ J dtdxcos{l^(j)i), (7) 

iGbosons 

which corresponds to the bulk Higgs terms. Canonical 
quantization of 5^ yields the Kac-Moody algebra 

[d,(j)iix),dy(j)j{y)] = i27:{K-')ijdJix ~ y). (8) 

For simplicity, when referring to a quasiparticle ai or 
its edge mode l^(f>i, hereafter we will most often simply 
specify only the charge vector 1. Besides, since we are 
mostly interested in the fate of the edge states, hereafter 
we shall refer to Eq. (7) or simply the cosine functions 
therein as our Higgs terms. And we shall from now on 
focus on the edge modes (f) exclusively. 



B. Edge Gapping Conditions in LRE Phases with 
Symmetry 

To gap out a bosonic edge mode, one needs to con- 
dense it at certain classical expectation value; however, 
the uncertainty principle due to Eq. (8) may prevent 
one from doing so. Any two bosons labeled by vectors l^ 



and Ih must satisfy the following canonical commutation 
relation, implied by Eq. (8), 

[lid^<j>i{x),lidycj,j{y)]^i27T{K-')jjlilidJ{x^y). (9) 

It is then clear that if boson can condense, the above 
commutator must vanish for a ^ b, i.e., l^K^^la ~ 0. 
Such a boson is called self-null. Furthermore, two 
bosons la and lb are only allowed to condense simultane- 
ously if they are both self-null, as well as mutual-null, 
namely Y^K~^li,=Q. 

We now summarise the necessary and sufficient condi- 
tion for a non-chiral LRE Abelian phase with symmetry 
to attain a symmetry preserving yet fully gapped edge 
state. 

Complete-Gapping Condition Given a non-chiral, 
Abelian, LRE phase characterized by a K matrix sat- 
isfying dim if = N G 2Z+ and \K\ > 1, with a global 
symmetry group Gg, in order that the edge modes in the 
phase can be completely gapped without breaking Gg, there 
must exist at least one complete set Bj of independent 
self and mutual null bosons, of which the Higgs terms are 
invariant under the action of Gg , namely, Vg £ Gs , 

laeBi 

^ ^ Ci„cos(;^p„(5)0-^ + aiJ (10) 

laGBi 

where each la G Bj appears at least once, ai^ is some 
arbitrary angle allowed by the symmetry transformation, 
whereas pij{g) is an abstract notation of the represen- 
tation of g € Gs. The completeness of the set follows 
from two criteria: first, any boson that is not in Bj is 
not mutual-null with at least one member in Bj; second, 
the set must consist of at least N/2 linearly independent 
charge vectors. 

If this is the case, all edge boson modes can be pinned 
at a vacuum expectation value without breaking the sym- 
metry, and we call the corresponding LRE phase with 
symmetry edge-trivial. 

This condition should be self explaining. It is adapted to 
the case of LRE phases with symmetry from its counter- 
part in Ref'^ for the SPT phases, where one finds detailed 
reasoning for the condition. Two remarks are in order, 
however. For a K matrix describing a chiral LRE phase, 
it would be impossible to gap all the edge modes, as there 
would always be excessive left or right moving bosons. 
Also, for each given K matrix model there may be more 
than one complete set Bj which by definition cannot be 
mutual null with each other, and any single set that is 
completely gapped is sufficient to gapping the edge. 
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C. Representations of the Symmetry 

Given any Lagrangian, one could look for the sym- 
metries that leave it invariant. In our case, the La- 
grangian comprises the Chern-Simons terms and the 
Higgs terms. However, we will work backwards in the 
program of studying non-trivial phases with symmetries. 
We would start with the K matrix theory with a fixed K 
matrix that has the correct degeneracy appropriate for 
the phase wc are interested in (a topological phase with 
\K\ = > 1), then exhaust all possible group action 
on the excitations for a given symmetry group Gg. An 
element of Gs acts on the anyons depending on the rep- 
resentations they fall into, as indicated in Eq. (10). We 
would consider a very general scenario where the group 
actions rotate a dyonic state in addition to attaching a 
phase to a charge or flux excitation, which can be imple- 
mented by the following pair 

p(5)-{M^^d<AS}, V.geG„ (11) 

where is some GL{N,Z) matrix and a constant 
A^-component vector, called a shift vector, such that 

K W'^KW = T]K, (12) 

where the second line follows from our requirement that 
K is fixed, and rj ~ 1 for unitary actions and —1 for anti- 
unitary actions such as time-reversal symmetry. But we 
shall not consider any time-reversal in this paper and 
hence set rj = 1 here onward. As reasoned in Secion 
II, this representation is in general a projective (or even 
more general non-linear) representation of Gg- In partic- 
ular, to be consistent with the fusion group of the anyons, 
the action of the identity e £ can transform individ- 
ual anyons but has to preserve the physical quasiparticles, 
namely the bosons, modulo 27r, i.e.,. 

l'B{{Wn'Uj + dcj,J) ^l's4>i (mod 2^), (13) 

where 1^ labels the charge vector of a boson. Then we a 
posteriori look for Higgs terms invariant under the above 
action, and finally check for the presence of any remain- 
ing ungapped edge modes, or a gapped edge that dynam- 
ically breaks the symmetry. 

One can immediately infer from Eq. (13) that the ac- 
tion of {W^,d(f)^} does not depend on G^; it is simply 
the set of transformations that is as exotic as possibly 
allowed by the fusion group. It is equivalent to a projec- 
tive representation of the identity of Gs- Recalling the 
discussion in Section II, when Gg is present, {W'^,d4>'^} 
furnishes a linear representation of an emergent gauge 
group Ng. Wc thus rc-emphasize here that a projective 
representation {W^ ,d(f)^\g G Gs} can be interpreted as 
a linear representation of a total group G that is an ex- 
tension of Gs by Ng, This will be played out explicitly in 
our examples. 

Here is how one computes {VF^,d^^}. Because Eq. 
(13) must hold for any boson, one readily sees that = 



1. Then, one needs to solve Eq. (13) for d^'' with = 
1 plugged in. One would obtain a vector with a number 
of integer parameters determined by the fusion group. 
Each choice of the parameters renders the corresponding 
d^'' a generator of the emergent gauge group Ng. Each 
such generator then determines what Ng is for this choice. 

With {W'^,d4>'^} in hand, one can then solve for 
{W^ ,d<p^} by solving a set of equations, which follow 
from the group compatibility conditions of Gg. This is 
a list of independent group multiplication relations that 
completely specifics the group structure. For a simple 
example, if Gs =^21 there is just one such compatibility 
condition, = e, £ Z2, which is then translated into 
the following equations: 

(T^f )2 = = 1 

W^iW^cl) + d03) + d(j)^ = Wcj) + dfj)" 
=^ M^^d^^' + d</)» = d0^ 

More complicated Gs may have more than one and more 
complicated group compatibility conditions. Among all 
solutions for , one can only choose those satisfying 
Eq. (12). The shift vector dcf)^ contains the parame- 
ters of d^'' and its own parameters in general. If the 
parameters in d^'^ are all switched off, it clearly implies 
that Gs stays unextended: G ~ Gs- Otherwise, one may 
need to work out the group compatibility conditions to 
determine precisely the type of group extension and the 
total group, for each choice of the parameters. Although 
Gs is an Abelian group in this paper, the actions by the 
projective representation of two distinct elements on an 
anyon do not necessarily commute. This is particularly 
true when at least one of the actions is represented by 
a nontrivial W^. This non-commutativity is ubiquitous 
in our examples, e.g., as explicitly discussed in Section 
VI B. 



1. Note of caution: residual gauge symmetry 

It is important to realize that the K matrix construc- 
tion suffers from a lot of redundancy. Different K matrix 
can describe identical phases, if they are related by rela- 
beling of the anyons, leading to, K X^KX, for some 
X € GL{N,Z). As aforementioned, we would be work- 
ing with specific K matrix which is known to describe the 
topological phases we are interested in. Even with a fixed 
K matrix, however, one is still haunted by the relabeling 
redundancy, because there are residual reparametrization 
X which keeps a given K matrix fixed. As a result, 
different global symmetry transformations may not be 
uniquely defined. For those related by X in fact describe 
precisely the same phase. More precisely two sets of 
transformations {W^%dcf)^^} and {W^%d(p^^} describe 
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the same physics if they are related by the following^ 
K = X^KX 

W"' = X-^W<''X, (14) 
d(j)'^' = X-\d(j) + A(l)^ T^9-A</.), 

A(/i is an arbitrary vector with integer components. 
These relations will help in locating the most convenient 
representative among equivalent solutions for and 



IV. SEP{Z2 X Z2,Z2) PHASES 

In this section we construct our first example of sym- 
metry enriched topological phases characterized by fusion 
group Z2 X Z2 and global symmetry Z2 . It is known that 
this fusion group is shared by two admittedly distinct 
models of topological order, the double semion model de- 
scribed by the K matrix ( -2 ) = '^'^z , and the Kitaev's 
toric code model defined by the K matrix (22)" '^'^x- 
We shall incorporate Z2 symmetry to these two models 
in order in the following two sub sections, then study the 
relations between these phases. 



A. SEP{Z2 X Z2,Z2) from the Semion Model 

As aforementioned, the double semion model is de- 
fined by the K matrix 2az, which is invariant under the 
GL(2,Z) transformation ±l2, and iciz- The quasiparti- 
cle content of this model is determined by the self statis- 
tics in (2) and described as follows in terms of the vectors 

Semions: 

The self statistics {K-^)ijlH'^ = ±1/2 (mod 2) 
demands h = h + (2m -I- 1), Zi, I2, m G Z. Hence, a 
generic semion is 

I5 = {l,l + 2m + 1), l^mf^'L. 

The elementary semions (of opposite chiralities) are 
thus SL = (1, 0)^ and sr = (0, 1)^. 

Bosons: 

The self statistics {K-^)ijfl^ = {If - ll)/2 = 

(mod 2) sets G 2Z. Thus, a generic boson 
takes the form 

lU = (2to, 2n), m,n £ 

which is an authentic boson because it has triv- 
ial mutual statistics with an arbitrary quasiparti- 
cle F = (^1,^2): 2{K-^)ijiyJ = 2{mk - nh) = 
(mod 2). The elementary bosons are (2,0) and 
(0,2). 



Bosonic bound states of semions: 

There are also quasiparticles consisting of both el- 
ementary semions that have bosonic self statistics 
but nontrivial mutual statistics with the semions. 
These take the following general form. 

iff, = (2m + l,2n+ 1), m, n e Z, 

where the subscript stands for bosonic bound 
states. 

Sets of independent condensable bosons: 

A condensable boson is a boson as defined above, 
with the additional requirement that its self statis- 
tics is identically zero, and not only zero modulo 
2tt. a condensable boson therefore has to satisfy 

{K~'^)i.j1^bIb = 2(m2 - n^) = 0, ^ m = ±n 

Multiple condensable bosons can condense at the 
same time only if their mutual statistics is also 
identically zero. In this case therefore, the 
two independent sets are given by {fc(2,2)} and 
{fc(2,-2)}, for aU fc e Z. 

Here we remark that fcrmions arc not in the quasi- 
particle spectrum of the double semion model because 
II — I2 7^ 2 (mod 4), V/i,Z2 S which disallows fermionic 
self statistics. Having listed out the quasi-particlc in the 
model, we can try to solve for all possible (projective) 
representations that are consistent with the fusion prop- 
erties of the quasi-particles when incorporating a global 
symmetry. The idea is that the identity element of any 
global symmetry must act trivially on each and every 
bosonic particle, a condition already described in Eq. 
(13). Yet this does not necessarily imply that it acts 
trivially on all quasi-particle excitations, although that 
is one obvious option. Therefore the first step we take 
is to solve for all possible non-trivial "identity transfor- 
mation" compatible with Eq. (13), which we label by 
{W'^,dr}. 

For the semion model, the solution to Eq. (13) is given 

by 

W^ = l, dcf)" = n{ni,n2f, ni,n2e{0,l}. (15) 

This means that we have altogether four different op- 
tions at our disposal. We can pick one to be the identity 
transformation for each independent symmetry group we 
introduce for each phase. Each such choice gives rise to 
an emergent Ng = Z2 when a global symmetry is incor- 
porated. As we will sec however, some of these difference 
choices could still potentially lead to the same phase. 
Next we have to solve for the rest of the symmetry trans- 
formations for a given symmetry group. For simplicity, 
we will consider incorporating a Z2 global symmetry here. 
This requires solving for the transformation correspond- 
ing to the single generator of the group, which we label 
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as {W^ ,d4>{g)}. Since = 1, this transformation must 
satisfy 

W<^dcj)3 + = d^'' (mod 2tt), 

{W^f ^^ = 1. ^^^^ 

There are several sets of {W^ ,dcj){g)} that satisfy the 
above. They are listed as follows: 

^±l,d0 = V2("i,n2)^, ni,n2G{0,l}} 
{W<^ = ±<j,,dct,^n/2{ni,n2f, ni,n2 G {0,1}} 

At first sight there are many possibilities. However, we 
note that when = — 1, the transformation dcp is not 
invariant under residual gauge transformation, and can 
be entirely gauged away. Therefore this choice corre- 
sponds to the same phase as the semion model without 
symmetry. Also we note that the choice for ti , t2 has no 
effect on the transformation of any bosons up to shifts 
of multiples of 2mT. Therefore it has no bearing on the 
allowed Higgs terms, and therefore different choices of 
which would lead only to the same phase. Therefore 
w.l.o.g we will consider the representative case where 
they are chosen to be zero. When is not the identity, 
the group action corresponds to swapping quasi-particles 
or anyons around. This will be considered in a later sec- 
tion. We will focus on = 1, and consider separate 
choices of ni,n2 in turn. 

Case Ia:S10: For later convenience, we label the phase 
for {71,1,77.2} = {1,0} by 510, where S stands for 
the semion model. The invariant Higgs terms are 
given by 

Sf' = cos(4m(0i - 02)). (18) 

As a result, the bosons with charge vector 2{2n + 
1)(1, —1) which is shifted by tt under dcf){g) would 
either acquire a vev due to the Higgs terms above 
and thus break the Z2 symmetry, or would remain 
gapless. 

Case Ib:S01: The case S'Ol works very similarly since 
it as far as the Higgs terms arc concerned, it is a 
relabelling of bosons by (pi — > —02, 02 — >■ — 0i- The 
edge therefore remains gapless. 

Case II: Sll: The invariant Higgs terms arc given by 
^1" - J2 ^™ cos(2m(0i - 02)) (19) 

Clearly, in this case all mutually condensable 
bosons in the set 2777,{1, —1} can all be simultane- 
ously gapped. Therefore the S'll phase has trivial 
edge. 



B. SEP{Z2 X 12,^2) from the Z2 Spin Liquid 

The K matrix taken as our starting point here is given 
by K = 2a X- Similar to the semion model we begin by 
listing all the quasi-particles: 

Self-commuting "boson": These are excitations that 
have bosonic self statistics; however, they do 
not have trivial mutual statistics with all other 
excitations, which is why they are labeled as 
"bosons" in quotes. Their charge vectors are 1"^ = 
{(771, n)|m, 77 £ Z, 777, = 7T, -|- 1 (mod2)}. Hence, the 
elementary "bosons" are (1,0) and (0,1). 

Fermionic bound states: 

These are the set of particles that have fermionic 
self statistics ~ tt (mod 27r). The charge vectors 
are given by 1"^ — \2n -f 1,2™ 4- 1}, 77, m G Z. 
The elementary "fermion" is given by charge vector 
(1,1). 

Bosons: These are true bosons with 2n7r mutual statis- 
tics with all quasi-particle exciations, and 2777r self 
statistics. Repeating precisely the same exercise as 
in the case of semions, we arrive at the set of charge 
vectors 

Is = (277,2777), 77,771 GZ. (20) 

Sets of independent condensable bosons: 

Straightforwardly, the two independent sets of mu- 
tually commuting condensable bosons are {(2r7, 0)} 
and {(0,277)}. 

Similarly to the semion model we can solve for all pos- 
sible "identity transformation". The solution is identical 
to that in Eq. (15). Consider again imposing a global 
Z2 symmetry on the Z2 gauge theory, we then solve for 
sets of transformation matrix {M^^, d</)(gi)}. The distinct 
solutions are 

{W^ = 1, d0^ = V2(ni, «2), m, 7^2 e {0, 1}} 

{W3 = ±a„d0.^ = V2(ni,n2), «i,7i2 e {0,1}}, 

(21) 

where we have already hidden a possible 7r(ti, <2) in d(^ 
under a rug for the same reason as before. Focusing on 
the solution on the top line — 1, the distinct phases 
are as follows: 

Case I: TIO: We adopt similar labelling of the distinct 
phases, and "T" is an allusion to the Toric code 
model due to Kitaev, which is a popular solvable 
model realizing the Z2 gauge theory. Here, the al- 
lowed Higgs terms are 

^i" = II cos(2?7702) (22) 

This clearly exhausts an entire set of mutually 
commuting condensable boson. Therefore the TIO 
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phase has trivial edge. The TOl phase is obtained 
by a relabeUing 0i O 02, and thus gives the same 
phase. 

Case II: Til: Here the allowed Higgs terms (from a sin- 
gle mutually commuting set of bosons) are 



cTll _ 
Of? — 



Cm cos(4m(/)i). 



(23) 



Bosons with charge vector 2(2n + 1)(1,0) there- 
fore either remain gapless or breaks the symmetry. 
Therefore Til has non-trivial edge. 



C. A (quasi)-group structure between the phases 



We consider superposing different phases with non- 
trivial edges found above in turn: 

Til X Til. Consider superposing Til phase with 
whose edge modes are denoted and another 

Til with edge fields The K matrix of the 

combined system is the direct sum of that of the 
constituent models. In this case therefore it is given by 
^TiiTii = 2(72; ® 2(7^ One allowed sets of Higgs terms 
within a chosen set of mutually condensable bosons are 
given by 



5i"xT"^^Cicos(2m(0i-0i))- 



E 



(24) 



Now we would like to discuss a (quasi)-group struc- 
ture that emerges by superposing the distinct phases with 
global Z2 symmetry we have obtained using the semion 
model and the Z2 gauge theory as the starting point. 

Our discussion of a group structure closely follows that 
in Ref'^. The basic idea there is that one can define a 
group product structure between two phases A and 
within a class of phases with a given symmetry, by stack- 
ing one on top of the other. The combined phase would 
generally allow for extra Higgs terms, gapping further 
edge modes. When a group structure is well defined, 
one could show that the combined phase, described by 
a new K matrix that is the direct sum of those of the 
component phases, can be transformed after appropriate 
reparametrizations, into a direct sum of a trivial SPT 
phase with a gapped edge, and another that is a memeber 
C of the original class of phases with the given symmetry. 
This allows one to identify a group product A® B = C 

There is a crucial difference between SPT phases and 
our LRE phases. In the case of SPT phases \K\ = 1, 
this is preserved as we superpose phases. This allows 
one to naturally dump the SPT phase whose gap is triv- 
ially gapped after we stack the phases. This is no longer 
the case when we have \K\ > 1. Therefore the group 
structure we are aiming for is not strictly a group. But 
consider the following situation: suppose we put two 
phases, A and B together, each with a non-trivial edge, 
and put them together exactly as in the procedure de- 
scribed above. Suppose also that there exists a relabeling 
of the bosons such that the new reparameterized K ma- 
trix becomes again a direct sum of two topological phases 
with the Higgs terms now diagonalized in each compo- 
nent phase, and that at least one of which has entirely 
gapped edges, and the other, called phase C, is recogniz- 
able as one of the phases we defined before superposing. 
Then there is indeed some notion of a group structure 
where the group product oi A x B = C, and that all 
phases with trivial edge are treated as the identity el- 
ement. As we will find below, such a group structure 
indeed exists, but the group product only closes if we are 
allowed to include both the S and K models in the group 
product. 



One can check that this exhausts the entire set of mutu- 
ally condensable bosons. The combined phase is left with 
a trivial edge. To display the group structure, we now 
considering a relabelling of modes given by the following 
conjugation K, X^ICX for some 5*^(4, Z) matrix X: 



X 



/I 
10-1 
10 10 

Vo 1 



(25) 



This matrix X leaves /Ctiitii invariant. However one 
can check that under this reparametrization where (f> 
'K~'^(j>, the group action = 7r/2(l, 1, 1, 1) after the 
transformation becomes 



X-^dcf) = 7r/2(l,2,0,l)^ 



(26) 



The entry with value 2 in the transformation vector above 
acts trivially on physical bosons, whose charge vectors 
consist only of components divisible by 2. Therefore, 
Til X Til is indeed the direct sum of two phases that 
we have already encountered previously: 



Til X Til = TIO® TOl - 1, 



(27) 



and each of TIO and TOl has trivial edge, and so we 
replace them by "1", the identity element. 

SIO X SIO. Following the same logic as before, by 
superposing two SIO phases, we arrive at the model 
A^siosio = 2(T2 © 2az- The allowed Higgs terms are 



■SlOxSlO 



J2 cos(2m((/.i - ct^],^ -cj^l- (/.!))+ 



mG2 



cos(2to( 



(28) 

Despite the appearance of two independent sets of Higgs 
terms, one can see that there are further mutually con- 
densable bosons that breaks the symmetry. They are 
{2m(0[-0)j)} and {2m{(j}l + . Therefore SIO x SIO 
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contains non-trivial edge. To make contact with the orig- 
inal 2 X 2 K matrices, we consider again a conjugation 
transformation of ICsiosio by X 



X 




(29) 



This transformation does not leave JCsiosio invariant. 
It is transformed upon conjugation into X^ICX = 
2a z © 2ax- Correspondingly the group action d(p^ = 
7r/2(l, 0, 1, 0) becomes 

X-^dcf) = 7r/2(3, -1, -1, 3)^ ~ 7r/2(l, 1, 1, 1)"^, (30) 

where wc use symbol ~ to mean that X^^drp, when act- 
ing on physical bosons is indistinguishable from the final 
transformation vector on the right. Therefore wc con- 
clude that 



SIO X SIO = Sll ® Til - Til 



(31) 



where Sll as we recall has trivial edge and we define our 
group structure that is only sensitive to the phase that 
has non-trivial edge states. 

SIO X SIO X Til ~ SIO X SIO X SIO X SIO. From 
the above, we can immediately conclude that 

SIO X SIO X Til ~ Til X Til - 1, (32) 

and that 

SIO X SIO X SIO X SIO - Til X Til - 1 (33) 

SIO X SOI. One can easily check that the combined 
phases allows the following set of Higgs terms 



jSlOxSOl 



^C^cos(2m(0i-4)) 



-C^cos(2m(0i, -</)!)) 



(34) 



which exhausts all mutually condensable bosons, and 
thus have a trivial edge. 

Therefore, one may be tempted to collect all the phases 
characterized by the fusion group Z2 x Z2 with Z2 sym- 
metry and arrange them according to the emergent group 
structure 



(35) 



SEP{Z2 X Z2,Z2) D Z4 = 

{1 - [TOO, no, 5*00, S'll], 5*01, Til - SOl'^, SIO}. 

Let us also clarify here that by taking the phases 
[TOO, no, 500, 511] to be the identity of the group struc- 
ture, it is not to be understood as identifying these 
phases. In fact as also emphasized in'^^ these phases can- 
not be connected smoothly without a phase transition or 
breaking the symmetry. We note that this collection of 
phases in Z4 do not include phases that involve non-local 
transformations of the anyons where = zka^ in the 
double semion model and = zt(Tx in the toric code 
model. As we shall see in Section VIA, there are ad- 
ditional phases whose edge always remains gapless, and 
that stacking them together never lead to a gapped edge. 



V. GENERALIZATION TO Zjv GAUGE 
THEORIES WITH GLOBAL Zm SYMMETRIES 

The discussion in the previous section over endowing 
the semion/ Z2 gauge theories with a Z2 symmetry can be 
readily generalized to the case of taking some (general- 
ized) Zat gauge theories and introducing Za/ symmetry. 

A (generahzed) Z^r gauge theory can be described by 
a K matrix of the following form: 



[N.l] 



N 
N 21 



(36) 



where N,l ^ Z, and ^ S {0, 1, • ■ • iV — 1}. They are in one- 
to-one correspondence with the Dijkgraaf-Witten lattice 
gauge theories, or equivalently the TQD models^^ 

Given K{N,l), one could readily obtain the general 
form of the physical bosonic excitations in the model. 
They are given by charge vectors 1^ = {li,h)'^ of the 
form 



li = Nnii , 



I2 — Nm2 + 2mi/, 



(37) 



where mi, m2 G Z. 

There are two independent sets of condensable bosons 



A:={m(0,iV)^}, B : {- 



Nr 



-{N,lf}r^,^. (38) 



by 



?cd(iV,0 

One could solve for the set of {W^ ,d(ff} . It is given 



1, 



27r 



2ln2 
N 



"-2 



(39) 



for Ui S {0,1,---7V — 1}. From Ref^ it is asserted that the 
allowed projective representations of the symmetry group 
consistent with the fusion algebra always take values in 
the fusion algebra itself. For any I, the fusion algebra is 
additively generated by d^'' with all possible choices of 
nil and 712 values, which can be straightforwardly derived 
as 



I = 



Zfcjv X Z^/fe, I ^ 



(40) 



where k N/gcd{2l,N) for N e 2Z, and k = 
N/ gcd(/, N) for N £ 2Z+ 1. On the other hand, since we 
are considering a single Gs = Z2 symmetry, the "gauge 
group" Ng that is involved in extending Gg is additively 
(mod 2tt) generated by a particular d^'' with a specific 
pair of ni and 77.2 values, as shown in the following equa- 
tion. 



N/xi 

yN, 
N/z, 



ln2 = 

2/an2 < N 
2lan2 > N 



(41) 
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where 

X = mm[gcd(ni,7V),gcd(n2,7V)], 
y^N/gcd{2ln2,N), 

z = mm[gcd(|rii - 2/r7,2/A^|, iV), gcd(n,2, iV)]. 

The corresponding transformation generated by the 
generator of a Zm global symmetry takes the form: 



= 1, d</>» = — 




(42) 



where here we are still focusing on cases that do not 
involve rotation of anyons, and iiG{0,l,---M — 1}. 

To determine whether the edge is gapped in each of 
these cases, we compute the transformation of bosons in 
each of the two complete condensable sets of bosons. If 
either set can be completely gapped, the edge is gapped, 
but otherwise remain gapless. 

The transformation of the bosons in each set is given 

by 



lie 



hj — 27rm 



27rr 



N^ti + Nlt2 + Nni - 712 
gcd(iV, l)M 



(43) 



The edge would be gapped if cither transformation 
vanishes modulo 27r with no further constraint on m or 
r. In general it would require specifying M and N and 
also the set {ti, rii} before one could determine if an edge 
has been gapped. Nevertheless let us illustrate in a few 
examples some representative cases. 



particularly simple forms: 
-) 

-{Nti+lt2 + ni-bi2/N), 



lie 



lie 



2'Km{t2 
27rr 



(45) 



5cd(iV,0 



and a non-trivial edge is formed if neither of the two 
sets of condensable bosons can be completely gapped out 
without breaking the global Zjv symmetry. From the 
behaviour of the bosons in set A, it is immediately clear 
that whenever 



772 = (mod N), 



(46) 



set A is gapped, independently of the value of ti and tii. 
In fact one can check that the value of ti is immaterial 
in the transformation of any physical bosons. There- 
fore they do not parameterize distinct phases and will be 
dropped from now on. 

Suppose that gcd{N, I) — x, so that we can write N = 
xa and I = xb for a, b relatively prime. Then the gapping 
of set B modes requires that 



77i — 6712 = (mod N). 



(47) 



Non-trivial edges therefore arise if Eqs. (46, 47) are not 
satisfied at the same time. This leaves, for each I a set 
of phases with non-trivial edges parameterized by 




(Si + 6n2)mod N 

n2 



Si, 772 e {i,---7V-i}. 

(48) 



C. A quasi-group structure 



A. M = 2, N = odd 

For simplicity, let us begin with a very specific exam- 
ple. In this case, we find that as soon as Nt2 + ^2 is 
even, set A is completely gapped, and thus the edge is 
gapped. Therefore we need only to consider what hap- 
pens if Nt2 + k2 is odd. In that case, we have to deter- 
mine if set B can be gapped. Let gcd(A^, /) =: x, such 
that N = xa and I = xb, with a, h relatively prime, the 
condition that all set B bosons are gapped is then given 
by 

a{Nti + ni) + b{Nt2 - n2) = Q (mod 2). (44) 

Recall that both a and Nt2 — 772 are assumed odd here. 
This suggests that if b is even, the edge is gapped when 
Nti -\- 771 is even, and for b odd, so should Nti -f 771. 

B. Special case: M = N 

In this special case, the above shift transformation act- 
ing on any one set of the condensable bosons take on 



For the general case where N > 2, we find ourselves in 
a large network of phases, and it is by no means obvious 
that the simple (quasi)- group structure that we find for 
N — 2, A/ — 2 that arises as we stack multiple phases 
on top of one another should also arise here. Rather 
than giving a complete survey of the matter, which seems 
much more complicated, we restrict our attention to the 
case N = M = 3 and I — and demonstrate, in this 
restricted scenario, that there is still a group structure 
existing between the phases. 

When I = 0, the phases with non-trivial edges are pa- 
rameterized by different d<f)^ as follows: 

d0S = ^(771,772)^, (49) 

where 771,772 G {0, 1,2}. We focus on phases with sym- 
metry here, as that with 771 = 772 = is equivalent to 
the usual topological phase without symmetry. Given 
that when 771 = (772 = 0), the corresponding phase can 
by fully gapped by condensing with the variable 3m(j)ji, 
we shall focus only on the phases labeled by nonzero 
Til and 712: there are, up to interchanging 771 and 772 
by renaming of quasi-particles three phases, given by 

(771,772) = (1,1), (2, 2), (1,2). 
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1. Stacking (11) with (12) or (22) with (12) 

In these two scenarios, we find that the edges of the ag- 
gregate phase can be completely gapped out. The Higgs 
term take the following form 

rneZ (50) 



The corresponding Higgs terms for 5*^ * takes the 
same form as the above, except that the signs of 
arc flipped. 



I and 



2. Stacking up three phases of the same kind 

It is not hard to check however, that stacking two 
phases of the same kind lead to non-trivial edges still. 
The next simplest option is to stack up three phases of 
the same kind. Consider for example stacking up three 
(11) phases. One can check that there is a complete set 
of Higgs terms that gap out the edge, given by 



S 



iieiieii 



E 

rni ,7712 ,m3 £2 



nil ,m2 ,m3 



xcos (^3mi{cbl + + 



(51) 



The same set of Higgs terms applies also to stacking three 
of the (22) phases or (12) phases. 



3. Stacking two (11) and a (22), or vice versa 

The above results already give us hints of a group 
structure. But to prove our point, we consider also this 
case. It turns out that this is again completely gapped. 
And one choice of the complete set of Higgs terms are 
given by 



c-lie22®ll 



mi,m2,m3i 



rni , 7712 , ma 



X COS (^3mi{(f>l + (1)1) 

+ 3m2(<^i + 0i) -f 3TO3(0fl, - 



(52) 



Having looked at the stacking above, we can recog- 
nize that the emerged group structure corresponds to a 
Z3, if we identity (11) and (22), which is justified from 
the fact that stacking three layers of (11), or two lay- 
ers of (11) with one layer of (22) both lead to a gapped 
phase. The (12) phase is the inverse of both (11) and 



(22), again pointing to identifying (11) and (22) in this 
group structure. (Wc note that we have not defined care- 
fully the procedure to preserve the ground state degener- 
acy \K\ as in the case of the Z2 gauge theory and doubled 
semion model. However, the similarity with the previous 
case makes it sufficiently evident that it should work very 
similarly here. ) 



VI. SEP PHASES INVOLVING THE 
ROTATIONS OF QUASI- PARTICLES. 

As already mentioned while we analysed the Z2 gauge 
theory and semion model in detail, there are interesting 
choices of symmetry transformation involving a transfor- 
mation matrix that is not the identity. Such a possi- 
bility was already explored in the K matrix construction 
of SPT phases without topological order"^. When there 
is topological order, such transformations have a partic- 
ularly vivid physical interpretation. 



A. A return to Z2 theories with Z2 symmetry 

Let us return to the Z2 gauge theory with K matrix 
2(7^ , and recall that the allowed choice of = cTx which 
implements a global Z2 symmetry on the theory. Its ac- 
tion on the (j> is accordingly (j) — ~^4>, which alterna- 
tively, acts on the charge vector 1 as 1 W^^l. Recall 
that l"'" = (10) corresponds to the "electric" excitation, 
and that 1-^ = (01) the "magnetic" excitation, this sug- 
gests that the action of = is precisely to exchange 
the anyons, implementing an electric-magnetic duality in 
this case. In fact, more generally, whenever ^ 1 it 
permutes the anyon excitation. Such a symmetry opera- 
tion is non-local and is not considered in Ref^. 

We note also that whenever ^ 1, only eigenvectors 
of could stay invariant, up to a sign (since deter- 
minant of = ±1). However, since is directly 
proportional to the K matrix (and its inverse) itself, it 
implies immediately that these eigenvectors cannot be 
self- null at the same time. In other words, no condens- 
able bosons could be left invariant by the Z2 gauge 
theory or the semion model. Therefore all of these phases 
have non-trivial edges, and no amount of stacking among 
these phases can lead to a gapped edge. In this case, 
can be either or —ax-, indicating a Z x Z group of 
phases with gapless states Note that we have relaxed our 
definition of a group structure here compared to our dis- 
cussion in Section IV C where a Z class is referring to the 
fact that we keep getting new phases as we stack phases 
on top without ever hitting a phase with trivial edge, 
although without defining a corresponding procedure to 
remove part of the system to preserve the torus ground 
state degeneracy |_?sr|, strictly speaking these extra phases 
may not belong to SEP {Z2 x Z2,Z2) . 

The same consideration applies equally to the semion 
model, except that an admissible choice of which 
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keeps its K matrix invariant is given by icTz, indicating 
also Z X Z group of symmetry enriched phases. 



B. More exotic examples: Z2 x Z2 symmetries in Z2 
gauge theories 

Such a global symmetry is considered also in Ref^. 
When the symmetry group is a direct product of groups, 
one could imagine that there are several relations among 
the groups. In the case of Z2 x Z2, it amounts to the 
following: 



9l 



1, .9k = 1, 
9l9r9l9r 1- 



(53) 



where gj^ and are respectively in the left and right Z2 
factors of the global symmetry. For each relation, one 
needs not have the same choice of {W,d<f)^} replacing 
the action of the identity, up to some consistency con- 
straints. Had we chosen, however W^^- = W^^ = 1 as 
in the previous sections, the only transformation has to 
come from the shifts dep. We would end up with the 
statement that the operators implementing gj^ and g^ 
necessarily commute. And the analysis that follows from 
taking W^^ = W^^i — 1 would be very much similar to 
what we have already considered in the previous sections, 
which we will not repeat here. 

The choice of W^^ = = 1 indeed does not ex- 

haust all the possibilities. Particularly, as we inspect 
the examples given in Ref^, models have been directly 
constructed where the symmetry transformations imple- 
menting gj^ and anti-commute. Before diving into a 
thorough comparison of the K matrix construction with 
other constructions, we would like to explore such a pos- 
sibility in K matrix construction, and specifically by un- 
derstanding the Z2 topological theories. 

Therefore, to construct a model such that the action of 
the generators gj^ and gj^ satisfy non-trivial commutation 
relations and at the same time allowing for the possibility 
that charges fractionalize, W^i- ^ W^R. Let us therefore 
consider K = 2ax and make the following choice 



(54) 



d0«-=7r(fi + ^,t2 + ^)^ 



d09« = ^(,5i,,52)^ 



where Ui^^^ correspond to the identity action dip'^ we 
choose for the group relation gl^^-j = 1, and that for 
consistency we require also that — n|^. 

One could now compare the action of (7^5^ and that of 
gj^gj^- They now lead to different shifts, which are given 



Si 



(mod 2), 



by 



d09i.9H = _(nf + 2ti + 2Si,n^ + 2t2 + 282)'^ 

I (56) 
d(p9a9L - -(n^ + 2t2 + 25i,n{ + 2ti + 282^ 

These relations demonstrate the following. First, that ti 
and t2 can now make a difference since they can deter- 
mine the eigenvalue of d0^^ under VF^^ = cr^. Second, 
it is clear that the action of g^ and gj^ on a fundamen- 
tal anyon (ie. (1, 0)"^ or (0, 1)"^) can be anticommuting if 
d0^^ is an eigen- vector of with eigenvalue —1. Never- 
theless, such a eommutativity is most natural when the 
representation of gj^ is in fact projective; otherwise, a lin- 



ear representation, which has 



0, would imply 



that the action of d<p^^^'^ on a fundamental anyon pro- 
duces a factor of exp(7r) while that of d<p^'^^'^ a factor of 
exp(— 7r), which are in fact identical. 

There is, however, one special situation where fraction- 
alization is not necessary for anti-commutative actions 
on a fundamental anyon, which is achieved by taking 

= 1, ^2 = and that 5i = 82 = \. In which case, 



^ 7r(2, 1)^ 

d^SflSi, ^ 7r(l,2)^ 



(57) 



One can see that each fundamental anyon acquires oppo- 
site sign under the action of g^g^ and gj^gi^- In this case 
also since iy^« is proportional to the K matrix and also 
its inverse, the edges cannot be trivially gapped. 



VII. SEP{Zi X Z4,Z2) PHASES 

We now consider the case where we incorporate a 
global Z2 symmetry into the topological phases described 
by the theories defined by a family of eight K matrices: 



/- 



K = 



V 




2 


Vo 



2 

2ni 


n2 





n2 
2 

2713/ 



The corresponding d(f)^'^ and d^^^ are then given by where 

(55) 



(58) 

ni, 712,713 e {0,1}. These K matrices all have 
l^l = 16, indicating that there are 16 quasiparticle types 
in theory defined by each such K matrix. If 712 = 0, it 
is clear that these K matrices turn out to be the direct 
sum of the 2 x 2 iiT matrices in Section IV; hence, we 
can infer that with Z2 global symmetry incorporated, the 
SET phases will be just those already found in SEP{'L2 x 
Z2 , Z2 ) . New phases with nontrivial boundary modes 
may thus appear only if rt2 is turned on, such that the 
K matrix is not block-diagonal. We arrange the three 
integers ni through ??3 into an array [711 71271-3] and use 
this to denote the eight cases to be studied. 
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A. Fusion and Gauge Groups 

111 this basis of the K matrices, a generic quasiparticle 
1 = ('i, ^2, ^3, ^4) has its components li and /a labeling the 
charges, while I2 and Z4 labeling the corresponding fluxes 
respectively®. The self statistics is 



1 



-hhn2^l\n^)^hl^ (mod 2), (59) 



which obviously can take values in {0,1/2,1,3/2}. 
Hence, there are 16 elementary quasiparticles all told, 
consistent with | det(if[oio])| = 16. Wc would not record 
all these elementary quasiparticles here but note that 
they can be obtained by allowing l\ through in Eq. 
(59) to be either or 1 and grouped by their self statis- 
tics. 

The fundamental quasiexcitations arc the two charges: 
ei = (1,0,0,0), 62 = (0,0,1,0), and the two fluxes: 
mi = (0,1,0,0) and ma = (0,0,0,1). These four fun- 
damental excitations all have the bosonic self statistics 
but not trivial mutual statistics with all other quasipar- 
ticles, as can be easily checked. But they can fuse to 
physical bosons. We would like to nail down the general 
charge vectors of bosons in terms of these fundamental 
excitations, which also allows us to read off the fusion 
algebra of the quasiparticles in this theory. 

Let \b = {Ii,l2,l3,h)'^ be a generic boson and 1' = 
(Z'l, ^3, ^4)"^ an arbitrary quasiparticle, their mutual 
statistics is 



+ + (^ + ^3 + kn^^, 



(60) 



which must be (mod 2). The terms in the above equa- 
tion are grouped as in the second row therein because 
the free variables are l'-^ through l'^, whereas li through 
are constrained such that the mutual statistics is trivial. 
Now that I'l through l'^ are free and independent, the four 
terms in the second row of Eq. (60) must be equal to 
(mod 2) individually. Wc then infer that the most gen- 
eral constraints on li through ^4 arc I2 ~ 26, ^3 = 2c—bn2, 
I4 = 2d, and li = 2a — dn2, where a,b,c,d € Z are free in- 
teger parameters. Quite naturally, these constraints are 
independent of rii and n^. Thus, the physical bosons of 
the theory take the following general form. 

1b = (2a - dn2,2b, 2c - bn2, 2d), a, b,c,de Z. (61) 



equivalent to the trivial particle = (0,0,0,0) in the 
fusion algebra, Eq. (62) leads to the following relations: 



ei X Gi = 0, 62 X 62 = 0, 



mi X mi X (62)"^ 



0, 



m2 X m2 X (ej 



(63) 



where the exponent is formal, meaning that (ej° = 
and (e^)^ = ei, i = 1, 2. It is straightforward to check 
that the fusion algebra ^?[„j,i2„3] of the 16 quasiparticles 
respecting the above relations turn out to be 



■S[nin2n3] 



Z4 X Z4 



7l9 = 1 



(Z2 X Z2) X (Z2 X Z2), ?12 = 



(64) 



These two fusion groups can also be verified by the 
projective representation {W'^,d4)'^} of the identity of 
whichever global symmetry to be incorporated, as we 
now show. Since this identity must preserve any boson 
up to a 27r shift, namely, l^^cj)! + Igdcji'j = 

(mod 27r), we immediately have = 14 and the fol- 
lowing constraint on d(ff 



--2ad(l)l+b{2d(l)l-n2d<j)l) 
+ c{2A(t)l-n2d(j)l) + 2dd<j3l 



(65) 
(mod 27r), 



where 1^ takes the general form in Eq. (61), and the 
terms are grouped by the free integer parameters a 
through d. As such, each term in Eq. (66) should be 
(mod 27r). Clearly, d(/)f and d(/)| can always be either 
or TT, and their value determines the possible values 
of d0| and dc/ij respectively. It is therefore not hard to 
write the allowed d(f)^ in a compact form as follows. 



1, d</,^ 



t2 + 



\ 



'-fn2 



U=u...A e {0,1}, 



(66) 

which readily generate additively (mod 27r) the fusion 
group Z4 X Z4 if 71,2 = 1 and the group (Z2 x Z2) x (Z2 x Z2) 
otherwise, as those in Eq. (72). Again for Gg generated 
by a single generator, the possible "gauge group" Ng in- 
volved in extending Gs is generated by a d<p'^ with one 
specific choice of ti through t^. There are only two pos- 
sibilities: 



Z4, tin2 ^ or t3?7.2 7^ 
Zo , otherwise 



(67) 



We can thus identify the following four elementary 
bosons: 

(2, 0, 0, 0), (0, 2, -n2, 0), (0, 0, 2, 0), (-712, 0, 0, 2). (62) 

The fusion algebra is generated by the fusion rules of 
the previously defined four fundamental quasiparticles, 
namely ei, 62, mi, m2. Since bosons are considered 



Case with [nin2n3] = [0n2 0] 



Seen from Eq. (66), 712 dictates whether the K ma- 



trix Ki 



nin2n3] 



has two decoupled blocks and thus the 



form of the fusion group. Since rii and play no role 
in the fusion group, let us set them zero, i.e., we have 
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[ni 7^2713] = [On20] = [712]. The K matrix in Eq. (58) 
becomes 



Ku 



(68) 

The G'i(4,Z) transformations that preserve A'jqxq] are 
the matrices as follows 








2 



















"-2 






"-2 


2 


0/ 



± 



I2 Qf521 

ad2i I2 
I2 

7^21 



7^21 
I2 



I2 (/3- 

/3<521 

(A -1)^21 

-12 



1)<521 

-I2 

I2 

A(52i 



(69) 



where I2J and ^21 are respectively the 2x2 identify ma- 
trix, and the 2x2 matrix (5 0)1 ^.nd a,/3,7, A G Z pa- 
rameterize an infinite family of these X matrices. 

According to Eqs. (61) and (66), we immediately see 
that in absence of global symmetry, all edge modes can 
be gapped by condensing either of the following sets of 
bosons, 



Ai 
A2 
A3 



= {(0,26,2c-n26,0)|fe,ce Z}, (70a) 
= {(2a,0,2c,0)|a,ce Z}, (70b) 
= {(2a-dn2,0,0,2d)|a,de Z}, (70c) 



or any set in the following two infinite one-parameter 



families of sets BpgQ and Bq 



B„ 



= |(-2cp,2fo, 2c-bn2,2bp) 

p £Q, b, c, bp e Z, 1^2-^ — cp e Z^, 
= |(-2(i(7,26,2feq-6n2,2d) 

gSQ, b,d,bqeZ, n2^-dqezy 



(70d) 



(70e) 



Note that if 712 = 0, the system is only a stack of two 
copies of the toric code model that is studied in our first 
example; hence, the above sets in Eq. (70) of indepen- 
dent, condensable bosons will recombine to merely four 
sets, each of which consists of one of the four combina- 
tions of the independent bosons respectively in the two 
toric code model with Z2 symmetry. 



1. Representations of the Z2 global symmetry 

Similar to previous examples, for Z2 global symme- 
try to be incorporated, we look for (projective) represen- 
tations {W^ ,d4>^}g^z^of the Z2 global symmetry group 
that transforms the fundamental fields but may allow cer- 
tain independent Higgs terms. We should first demand 
that for all 5 £ Z2, (W^)^ = I4 and {Ws)^ K[„^_]W3 = 
K. The latter condition guides us to find the correct 



matrices from the X matrices in Eq. (69); hence, 
we obtain = ±14, Xp, X^. We are interested in 
inequivalent transformations, and since = X^ 
and = X^ are are related by a G'L(4,Z) transfor- 
mation preserving the K matrix, as X^^X^X\ — —Xjj, 
they arc in fact equivalent and will not be considered 
separately. Moreover, for any value of 7, one can always 
apply a GL(4,Z) transformation by certain X matrix in 
Eq. (69) that preserves the K matrix, while keeping the 
form of d<p^ in (66) up to redefinition of the parameters 
ti through ti, to set 7 = in X^. Thus, we conclude 
with the inequivalent transformations 

= ±14, = X^=o = ±l2 ® (T„ (71) 

where ax is the usual Pauli matrix and CS) the usual matrix 
tensor product. Note that the matrices with a -|- sign 
and a — sign in the front are not equivalent to each other 
under the transformation in Eq. (14). 

Before we proceed to nail down the corresponding d(f>^ , 
let us remark on the behavior of = ±l2 ® <Tx- The 
action of on a quasiparticle 1 is given by {W^)'^l, and 
to manifest the physics we let 1 = (ei, toi, 62, ^2), where 
ei and mi (e2 and 1112) are respectively the charge and 
flux associated with the first Z2 (second Z2) gauge group 
of the total Z^ x Z2 gauge group. Then for = I2 O Ua; 
we have 



{Wsfl^{±l2®axfl = ± 



(e2\ 

7712 

ei 
\mi/ 



(72) 



which signifies a non-local exchange of the two types of 
dyons, (ei,mi,0, 0) and (0, 0, 62, "12) respectively of the 
two Z2 sectors of the gauge group. Such a non-local 
exchange transformation by the global symmetry is ev- 
idently beyond the scope of symmetry fractionalization, 
as also reported in Ref^. Note that this exchange trans- 
formation exists for any choice of [nin2n^] , even if n2 = 0. 

We now solve for d</>^. Since in any extension of Z2 by 
Z2 X Z2, the latter exists as a normal subgroup; hence, 
the group compatibility conditions demands that 

(1 + WB)A(i)^ = dcf)'', (73) 

for any dcf)^ in Eq. (66). We solve the above equation 
for d<p^ in different cases of . 

(i) = — 14. Equation (73) has the unique, in- 
equivalent solution dcf)^ = and ti ^ t2 = <3 = t4 = 
must be set in dcp'^. Since cos(Z^d0|) is invariant under 

= —1, the global symmetry Z2 does not transform 
the quasiparticles at all, implying that the edge modes 
can be completely gapped out, resulting in a boundary- 
trivial phase that is identical with the phase without the 
global symmetry. 

(ii) = I4. The solution of Equation (73) clearly is 



d0f = TT 




-d0^ p,^i,...,4 e {0,1}. 



(74) 
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This nontrivial shift vector in general prevents the edge 
modes from being fully gapped, as it forbids any of the 
sets of independent variables in Eq. (70). Special cases 
do exist, e.g., ti = = Q in dcj)^ would allow the entire 
set A2 to condense, resulting in a boimdary-trivial phase. 
Nevertheless, a thorough study of all boundary-nontrivial 
phases and their quasi-group structure in this case turns 
out to be rather complicated because wc lack of a con- 
venient and systematic algorithm for computing the new 
sets of independent, condensable bosons in a stacking of 
many phases for large-size K matrices. While we are not 
able to unveil the full quasi-group structure of the phases 
in this case, we do have a partial result to summarize as 
follows but fill the details in Appendix A. 

A study of how the independent bosons in Eq. (70) 
transform by the shift vector in Eq. (74) show that the 
relevant parameters in Eq. (74) are ti, T2 = t2 — '^2^3, ^3, 
and T4 = ti — n2Pi , where new parameters T2 and r4 arc 
defined in terms of the old ones. As such, our experi- 
ence tells us that we can label all possible phases by the 
values of the string [tiT2t3T4], leading to 16 phases. Tab- 
ulated in Appendix A, 12 out of these 16 phases actually 
have fully gapped edge state without symmetry break- 
ing. There are four edge-nontrivial phases remaining in 
Eq. (A4) with non-trivial edge: 

[tlT2t3T4] =[1010] 

[1011] 
[1110] 

[1111]. (75) 

We have not explored the quasi-group relations between 
these four phases, which gets cumbersome as larger K 
matrices are involved. This should be worth a future 
attempt. 

(iii) = ±l2 (8) da;. In this case, one can apply 
the equivalence transformation in Eq. (14) first to turn 
arbitrary into a common simpler form by remov- 

ing any redundancy. It is not hard to show that by 
choosing X = 1 in Eq. (14), for any dcj)^ , one can 
always find a shift to eliminate the first two com- 
ponent of the dcf)^ , without affecting . As such, one 
can assume that in general, dcf)^ = (0,0, x,y)'^ , where 
X and y are to be solved. The equation above now be- 
comes (±a;, ±y, x, y)'^ = d<p'^ ^ which is soluble only when 
ti = zbia in dcj)^ , leading to 

d0« = iTT , (76) 

\t2 + 1^2/ 

with constraints ti — its and i4 = ±t2 on d^'' enforced. 

Interestingly, however, since this does not yield 
any nontrivial shift to the boson variables in the set A2 
in Eq. (70b), as /^^dc^if = ±27rc<i = (mod 27r), VU^ e 
A2, one can gap out all the edge modes by condensing the 
independent Higgs terms constructed from the bosons in 



set A2 as follows. 



^Ca,c / dMx{ cos[2(a0i -I- C02)] + cos2(a(?!)2 c0i)]}, 

(77) 

where 0^2) ^'^^ respectively the electric edge modes asso- 
ciated with respectively the left and the first and the sec- 
ond Z2 factors of the Z2 x Z2 gauge group. Therefore, de- 
spite a nontrivial exchange of the two quasiparticle types 
under = ±l2 ® (Jx and even symmetry fractionaliza- 
tion due to the nontrivial d(/)^, the corresponding phase 
remains boundary-trivial. 



VIII. BEYOND CENTRAL EXTENSION 

Our examples in the previous two sections demonstrate 
some novel features in the transformation properties of 
anyons when one relaxes a crucial requirement imposed in 
Rcf^, namely that the symmetry operators transforming 
the anyons have to be local. In the previous two sections, 
we have provided several examples where the exchange 
of anyons, a glaringly non-local transformation, can give 
rise to more exotic phases, some of which for example has 
already been reported in Ref^. 

There is another important class of phases which also 
generally involve non-local transformation of the anyons. 
Reiterating Rcf^ , restrictions to local transformations has 
led to a classification of phases via different allowed pro- 
jective representations consistent with the fusion rules. 
That, in other words, is equivalent to a classification of 
different central extensions of the global symmetry group 
Gs by an Abelian gauge group Gg that is taken to be the 
fusion algebra Gg = 5^ of the topological phase on top 
of which global symmetry is built^. The restriction to 
central extensions has been raised to include more gen- 
eral group extensions*. In this case, the global symmetry 
group becomes the quotient group Gs = G/Gg, and the 
gauge group is the normal subgroup of a bigger group 
G. Different phases correspond to different choice of the 
total group G for given Gg and Gg . In this construction, 
Gg is no longer the center of the group G, and so one does 
not expect the group action of Gg and Gg to commute. 
From the previous sections therefore, it almost immedi- 
ately follows that such group actions necessarily involve 
exchange of anyons. In fact, the examples in the previ- 
ous sections can be understood within this framework of 
general group extension. 

In this section we would like to make such a construc- 
tion within the K matrix framework more explicit, and 
illustrate these principles using a particular set of exam- 
ples, where the total group G is chosen to be one of the 
dihedral groups for some odd N . 
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A. Step 1: Obtaining a linear representation of G 

The virtue of identifying a total group G in the classifi- 
cation of phases with symmetry is that any non-trivial or 
non-linear transformation under the action of the global 
symmetry group Gs can be reduced to a simple linear 
group action in a suitable G. Here, wc will focus our at- 
tention on realizing G = for N odd via A'- matrices. 

In Dn one can specify each group element by a pair 
{A, a), where A = ±, and a = {0, 1, ■ • ■ , — 1}. Group 
multiplication between two such pairs is given by 

{A, a) ■ (B, b) = (A.B, {Ab + a)„,odJv). (78) 

A representation of each group element {A, a) is given 
by (VF(^',d0"), where = 1 and W'-~^ = M^^"' 

for some non-trivial VF^"' that keep the K- matrix in- 
variant, and that (VF'^')^ = 1. On the other hand, 
d(/)° = ^dcf)^~\ where d^*^^-* is ancigcn-vectorof 
with eigenvalue —1. The aggregate action of {A, a) ■ {B, b) 
is then given by 

W^^\W^^U + d(f>'') + del,'' 
= M^(^)M^(s)^ + VP^(^)d</,'' + d0". ^^^^ 

Since d</)^ is an eigenvector of W^^^ with eigenvalue A, 
and that components of (j) are defined only up to multiples 
of 27r, we conclude that the above representation is a 
faithful representation of Djv- 

In the special case where K = Nox for example, W^^^ 
can be chosen to be W^~'^ = Ux, and d^*-"-* = (1, —1). 

B. Step 2: Identifying the normal subgroup with 

the gauge group 

Having constructed a linear representation, we would 
then have to identify a normal subgroup Gg of the total 
group G such that the group action of Gg is taken to 
be unphysical. In other words, Gg is taken as some kind 
of gauge group that does not have any visible effect on 
any physical, or gauge invariant, excitations. Therefore, 
admissible Gg is strongly restricted by the fusion group 
^. In fact they are embedded inside ^. In other words, 
the normal subgroup Gg can only be chosen whose group 
action on physical bosons in a K matrix model is trivial, 
i.e. 

l5(WV + dct)3) = ll(t> (mod 2tt). (80) 

Now returning to our dihedral group D^. Suppose we 
would like to pick the Z^v normal subgroup as our gauge 
group. This Zn subgroup consist of pairs {A — +,a), 
where the first component is The group action is 
then given by {W^~^^ = l,d^''}. This can be admitted 
as a gauge group only if l|Jd0° = (mod 27r). This 
already suggests that Zn must be at least a subgroup of 
the fusion group. For example K = Nax, where = 



N{ni,n2), for any S Z, and dcf)" = 27ra/iV(l, -1)^ , 
indeed we have l^Jd^" = (mod 2tt), and therefore we 
are allowed to take Zjv to be a gauge group. 



C. Step 3: Implementing the global symmetry 
group 

The global symmetry Gg = G/Gg. In this case where 
G = Dn and Gg = Zjv, Gg = Z2. The group elements 
of Gs are the cosets of G w.r.t Gg. The identity element 
of Gs is the coset which is in fact spanned by the nomral 
subgroup Gg itself. In this case therefore, it comprises 
all the pairs (-l-,a). Other cosets are generated by the 
normal subgroup by left multiplication gg ^ g, gg £ Gg, 
and g S G. We note that right multiplication would yield 
identical cosets for normal subgroups Gg . The other non- 
trivial coset corresponding to the non-trivial element in 
Gs = Z2 is the set of pairs (— , a). 

Now, the final step is to pick any representative in 
one of the non-trivial coset, whose group action is now 
interpreted as that of the global symmetry group. It 
automatically acts non-linearly on the anyons. Its action 
is closed as a group, up to group action of Gg which is 
now so aligned with the fusion algebra ^ that physical 
bosons transform trivially. For our example at hand, we 
can take the generator of Gs = Z2 for K = Na^ to be 
{cjx,^^{l,-lf}, for any a G {0,1,- --^ - 1}. Note 
that the shift 2'Ka/N{\, —1) on any boson lead to shifts 
proportional to 27r. Therefore we need only to worry 
about W'^^\ Given that = Ux, it immediately 

reduces to a situation we have encountered already in 
subsection VIA, where not a single edge mode can be 
gapped as we continue stacking; hence, in this case, the 
quasi-group of the phases is Z. 

We note that the idea of central extensions work in 
precisely the same way, except that the group action is 
restricted to be commutative. Here we demonstrate how 
a non-trivial group extension can be implemented within 
the framework of K matrix construction. 



IX. COMPARISON WITH OTHER WORKS 

Endowing topological phases with symmetry is a novel 
and important question that has been a subject of much 
interest recently. In the previous sections, we have pro- 
vided yet another construction of these phases based on 
K matrix. Among the scenarios we have studied, various 
have already been discussed in the literature. We would 
therefore like to make a comparison with known results. 



A. Comparison with Ref^ 

To begin with, we comment on the relationship of our 
work with that in Ref^. In Ref^, the main targets are 
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Abclian topological phases endowed with global sym- 
metries whose action is localized near the vicinity of 
the anyons excited in the system. In those cases, it is 
demonstrated that the anyons can transform under pro- 
jective representations of the symmetry group concerned. 
These projective representations are consistent with fu- 
sion rules: namely that the identity element must act 
trivially on any physical bosonic excitations, even if the 
bosonic state is a composite of fused anyons. This con- 
strains the possible projective representations allowed for 
individual anyons. There are limited choices of how the 
identity element of the symmetry group can act on any 
anyon. In our explicit construction via K matrices, it is 
clear that the requirements on d^*^ coincide with the dis- 
cussion of allowed action of the identity element. Among 
all the specific cases we studied, of ^ = Z^r gauge theo- 
ries and their twisted versions such as the double scmion 
model for N — 2, every single consistent choice of projec- 
tive representations as dictated in Ref^, which are classi- 
fied by H^{Gs,^) is realized in our constructions. 



B. A comparison with Ref* 

It is also of interest to compare our work with Rcf*. 
In Ref^ it is proposed that a systematic construction of 
topological phases with symmetry is to consider topo- 
logical terms of SPT theories with symmetry group G, 
whose normal subgroup N is subsequently gauged. Such 
a theory should describe a topological phase with global 
symmetry given by the quotient group Gg = G/N. Spe- 
cific examples where G = Z4 and separately Z2 x Z2 arc 
considered, in which N and Gg are both Z2 in each case. 
Therefore, these two possible G's correspond to differ- 
ent (central) extensions of the global symmetry group. 
Moreover, for each choice of G there are several choices 
of topological terms, classified by H^{G, U{\)). They led 
to many different possible phases. One distinguishing 
feature between these different phases constructed is the 
braiding of excitations around magnetic charges of the 
global symmetry group ~ by magnetic charges, they re- 
ally correspond to multi-valued field configurations with 
branch cuts that end at a branch point. In the K ma- 
trix construction however, all field configuration is single- 
valued, and these extra braiding statistics are invisible to 
us. If we ignore them, then there is a one-to-one corre- 
spondence between the phases we constructed and the 
phases studied there. In the case where G = Z4, there 
are four phases constructed in Ref^, which is parameter- 
ized by a topological term with coefficient m which can 
take values in {0,1,2,3}. The correspondence with our 
construction is as follows: 



m = 

m = 2 

m = 1 

m = 3 



KIO 
Kll 
SOI 

sm 



On the other hand, when G = Z2 x Z2 there are eight 
phases with three independent topological terms, each 
parameterized by a coefficient = {0, 1}, for i ~ 
{1,2,3}. It is demonstrated there^ that fractionalization 
occurs if and only if 712 is non- vanishing. There arc four 
phases therefore that admit fractionalization, and again 
each of them directly corresponds to our construction: 



(010) 
(110) 
(Oil) 

(111) 



A'Ol 
A'Ol 
S-ll 
S-ll 



(82) 



(81) 



We note that in the above, two phases arc mapped to 
the same K matrix phase because as emphasized already 
these phases differ only if magnetic charges of the global 
symmetry is visible, which they are not in our K matrix 
construction. It is perhaps also surprising that here all 
the phases have trivial edge excitation! 

We finally note that the general proposal of Rcf* gen- 
eralizes the classification of phases via central extension 
to a general group extension of the global symmetry 
group Gs, although explicit examples considered there 
lie within the framework of central extension. In this 
work we provide first such examples of a general group 
extension realizing the proposal in Ref^ . 



C. Comparison with Ref'^ 

Finally, we would like to comment on the relationship 
of our work with that of Ref^. In Ref^, a (lattice) topo- 
logical gauge theory as defined in Ref^^'^* is taken as a 
starting point, whose gauge group G is chosen to be a 
direct product of Q and Gg ■ The action amplitude of the 
theory is characterized by difference choice of "topolog- 
ical terms" i' which are group cocycles in H^{G,U{1)). 
The G,, part is then ungauged, by restricting field con- 
figurations to be pure gauge, i.e. For each field config- 
uration where degrees of freedom sit on the links of the 
lattice, each of which labeled by a pair {gg,gs), where 
Qg & Q and gs G Gs, each gs at a particular link in the 
collection of degrees of freedom can always be written as 
gs = SisJ^ , where Sij S Gs and i,j label the vertices 
connected by the link concerned. In otherwords, each set 
of link variables {gs} can be replaced by a set of vertex 
variables {si}. Q however remains gauged, supplying the 
long-range entanglement needed in a topological phase. 
It is observed in various explicit examples that the pure 
electric excitations always transform linearly under Gs, 
and that by picking different u magnetic or dyonic ex- 
citations of Q can transform non-linearly. In particular, 
anyons can transform in different non-trivial represen- 
tations of the global symmetry group Gg, including the 
projective representations as discussed in Ref*, but not 
restricted to them. 

We specifically wish to comment on two of our exam- 
ples which are motivated by observations in Ref^. Be- 
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fore that, one should note the role played by Q in these 
topological gauge theory constructions. In particular one 
should be cautious and observe that the residual gauge 
group Q is not to be confused with the fusion algebra 
^ between all the anyons. As already emphasized in our 
overview of the paper, Q is the "deconfined" gauge group, 
and it is (subgroups of) ^ that is often being identified 
as the gauge group in Rcf^, which we have denoted Ng 
throughout most of our paper. It is however as expected 
that in these models, Q forms the sub-fusion algebra in- 
volving only the pure electric excitations. Without going 
into technical details, we can identify the K matrix de- 
scription that corresponds to Q = 1^2 and Gg = Z2 x Z2, 
which is discussed in section VI B. We found examples 
where the action of the generators of the two Z2 symme- 
tries anticommutc. In section VII, we also found exam- 
ples corresponding to = Z2 x Z2 and G,, = Z2, where 
the group action of Gs is manifestly non-local and ex- 
change the gauge electric and magnetic charges between 
the two Z2 gauge groups. Both of these cases are con- 
sidered in Ref'' and where such novel transformations 
are also observed. As already mentioned above, while 
these novel transformations are absent among pure elec- 
tric charges in Ref*", they are ubiquitous in the K ma- 
trix construction, which naturally provides the flexibility 
to incorporate non-linear transformations on any excita- 
tions, as long as they are consistent with the full fusion 
algebra. For this reason, our specific choice of the group 
action on the electric and magnetic charges may not gen- 
erally coincide with the examples in Ref^. We believe 
that the distinction between electric and magnetic exci- 
tations there is an artifact of the topological gauge the- 
ory construction. Given the direct connection between Q 
and the fusion algebra of purely electric excitations, one 
realizes that to achieve non-trivial transformations also 
among pure electric charges in the framework set forth 
in Rcf^, one is compelled to consider topological gauge 
theories where G is taken to be a non-trivial extension 
of Gs by Q other than the direct product that is being 
considered^. 



X. DISCUSSIONS AND OUTLOOK 

In this paper, we have been studying intrinsically topo- 
logical phases endowed with a global symmetry - the 
symmetry enriched phases as we dubbed - aiming at 
their classification and edge-state properties, by means 
of the celebrated K matrix formulation of effective the- 
ories of Abelian topological order. While methodologi- 
cally we extends the application in Rcf'^'^ of K matrices 
in classifying SPT phases to LRE phases with symmetry, 
we systematically adapt and integrate several principles 
imposed in Ref"^'^"^'^^ , particularly of how a global sym- 
metry may transform the anyons in a topological phase. 
These principles and the K matrix method guide us 
to constructing examples of symmetry enriched phases, 
along with clarifying a few important conceptual ques- 



tions, particularly the roles of various different groups 
play in classifying different phases. As noted in Ref^ it 
is the fusion group ^ of the anyons under consideration 
that constrains the action of a global symmetry Gs, in 
a way such that the identity of Gs acts trivially on any 
physical bosons although it may transform any individual 
anyon exotically, in which case the anyons may undergo 
local symmetry charge fractionalization and perhaps ac- 
companied by non-local transformations such as anynon 
exchange. The "gauge group" Ng involved in extending 
Gs is a subgroup of ^ that is projected (as the kernel 
of the projection map) into the identity of Gs and thus 
preserves the bosons, which indicates the existence of a 
larger group G that contains Ng as its normal subgroup 
and Gs is its quotient group G/Ng = Gs- Therefore, two 
different Gs actions do not commute in general, nor does 
the action of Gg and that of Ng, as shown in some of our 
examples. 

The K matrix approach offers a convenient way of ana- 
lyzing the relations among the symmetry enriched phases 
by stacking the phases, in the sense of arranging the K 
matrices respectively characterizing the phases into a di- 
rect sum and the corresponding Gs representations in a 
direct sum in the same order. In the examples we have 
shown, the various symmetry enriched phases for a given 
^ and Gs constitutes a quasi-group structure. In partic- 
ular in the case with 5^ = Z2 x Z2 and Gs = Z2, as ex- 
plained in Section IX B, the phases in the corresponding 
quasi-group are identified with the phases under the same 
setting in Ref* . To emphasize the prominent role ^ plays 
in these symmetry enriched phases we label the phases 
accordingly as SEP{'S,Gs)- This notion not only covers 
the symmetry enriched LRE phases but also embraces 
the SPT phases: the fermionic SPT phases with a given 
symmetry Gs comprise SEP(7j2tGs) because the fusion 
group of fcrmions is Z2, whereas the bosonic SPT phases 
all fall into SEP{'Li, Gs) because bosons have trivial fu- 
sion group Zi . 

Most of the examples we constructed are inspired by 
Ref'', but there are important differences that should be 
noted in our construction and discussion. First, we have 
carefully defined the notion of "gauge group". In par- 
ticular, similar to Ref^, it is what we denoted Ng that 
is pertinent: whereas in Ref^, the term "gauge group" 
refers exclusively to what we have denoted Q Second, the 
constructions in Ref'" gives rise only to flux fractional- 
ization; however, the K matrix method treats charge and 
flux on an equal footing, naturally allowing charge, flux, 
and dyons to fractionalize simultaneously. 

The K matrix method has another virtue: it enables us 
to study the fate of the edge modes explicitly, obtaining 
the condition when a phase may have gapless edge modes 
protected by the symmetry. Seen in the examples we 
constructed, symmetry charge fractionalization or more 
exotic transformations of the anyons under global sym- 
metry in a LRE phase is neither a sufficient nor a neces- 
sary condition for the phase to possess non-trivial edge 
states. Although we do not know if these phases that 
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have trivial edges, despite displaying exotic transforma- 
tions under the action of the global symmetry, may still 
be adiabatically connected to an LRE phase without any 
symmetry, as far as the edge property is concerned, in 
the quasi-group of all phases in a given SEP{^, Gs), we 
may treat those phases having a trivial edge on an equal 
footing with the phase with the same fusion group but 
without the symmetry, as if they arc projected into the 
identity of the quasigroup. 

Our first example, i.e., SEP{7^2 x ^2, ^2), is also partly 
discussed in Ref^, which already exemplifies the impor- 
tant role the fusion group of anyons plays. We realized 
every phase that appears in the classification in Ref^. 
Our construction however also involves phases that do 
not appear in^, when we allow for non-local group ac- 
tions that exchange anyons. Furthermore, in Ref^, the 
classification of symmetry enriched phases is equivalent 
to the classification of the central extensions of Gg by 
Ng] however, our examples also include a non-central ex- 
tension of Gg = by Z2 to the dihedral group D3, as 
anticipated in Rcf^. 

Inspired by^^, having observed that the various groups 
involved in characterizing a symmetry enriched phase are 
related hy ^ D G D Ng and G/Ng = Gs, and that Ng acts 
trivially but Gs acts nontrivially on the condensed edge 
modes, we are encouraged to redraw our picture of sym- 
metry enriched phases as an example of the Hopf symme- 
try breaking, first proposed and phrased in Ref^*'^^'^^ to 
account for anyon condensations, generalizing Landau's 
symmetry breaking. This new paradigm of generalized 
symmetry breaking may become most suitable to cope 
with the non-Abelian anyons endowed with a symmetry. 
We shall report our detailed studies clscwhcrc^^. 

Let us close with a discussion of interesting questions 
that should be more thoroughly addressed in the future. 
We now describe them briefly below. 

1. While we have a detailed analysis of SEP{Z2 x 
Z2, Z2) that probably exhausts all the phases in the class, 
our treatment of other examples requires further analysis. 
Particularly it would be of interest to explore whether a 
qusi-group structure can be generally defined. At the mo- 
ment it appears that the order of any such quasi-group 
in Zjv gauge theories grows at least linearly in iV, which 
makes an analysis very cumbersome quickly. A more effi- 
cient method is necessary for a thorough understanding. 

2. As far as clarifying a group structure of Zjy gauge 
theories are concerned, there is another specific question 
to be addressed. In our first example, we have seen two 
different models with topological order, i.e., the double 
scniion model and the toric code model, which share the 
same fusion group. When the same symmetry group is 
incorporated, they together lead to a set of symmetry 
enriched phases belonging to the same quasi-group: in 
particular the Z2 symmetric double semion model acts 
as a generator. These two models are actually described 
by the set of K matrices ( 2 2'ri ) with a single paramter 
ri = ±1, which defines the double semion model when 
n = I and the toric code model when n ~ 0. In more 



general cases, a class of different models of topological 
order can be specified by a multi-parameter K matrix. 
For instance, the K matrix in Eq. (36) characterizes 
respectively N models respectively for the N values of 
the parameter I = 0, . . . , N — 1. These models have dif- 
ferent fusion groups as in Eq. (40). Ultimately, this 
parameter I labels the N 3-cocycles in the cohomology 
group _ff'^(Z„, J7(l)) that classifies the corresponding N 
models, we are not able to answer at this moment the 
question whether the symmetry enriched phases charac- 
terized by respectively the fusion groups in Eq. (40) with 
the same symmetry group Gs would belong to the same 
quasi-group in a nontrivial way, as opposed to simple 
direct product of the quasi-groups characterized respec- 
tively by the N fusion groups and Gg. We arc not able 
to answer this question in general either and hence leave 
it for future exploration. 

3. We have considered only discrete gauge groups and 
unitary symmetries in this paper. It is of interest to 
construct more cases with continuous symmetry groups 
Gs and also those involving time reversal. 

4. Having observed non-local transformations of quasi- 
particles under Gs, e.g., the dyon exchange discussed be- 
low Eq. (72), and since non-locality is rather intrinsic to 
non-Abelian anyons, we look forward to extending our 
studies to the interplay between non-Abelian topological 
order and global symmetry. Unfortunately, this is be- 
yond the reach of the K matrix formalism and thus begs 
for new approaches. 

5. A recent paper by Vishwarnath and ScnthiP" found 
that some symmetry enriched topological phases in 2 + 
1 dimensions can only exist as the boundary of some 
SPT phase in 3 + 1 dimensions. Wc have realized some 
new phases based on general group extensions using the 
K matrix. It would be interesting to understand if the 
K matrix or strictly 2-1-1 models can exhaust all the 
phases based purely on consideration of group extensions, 
or whether some extra phases are again only realizable 
as boundaries of higher dimensional non-trivial phases. 

6. In the last stage of preparing this manuscript, we 
noticed a very recent paper by Levin'^^ that studies the 
conditions that allow for gapless edge states in a pure 
Abelian, non-chiral topological order without any global 
symmetry. It turns out that non-trivial edges can appear 
and that they are protected by the quasiparticle braid- 
ing statistics in the bulk, instead of by any symmetry. 
One such example is the v = 2/3 fractional Quanthum 
Hall system. The topological phases studied in our pa- 
per however have fully gapped edges in the absence of 
symmetry. It is of interest to extend our investigation 
to incorporating global symmetry in these novel phases 
discussed in Ref'^^. 

As we finish our paper, we were brought to the atten- 
tion of the work of Lu and Vishwanath'^^ which contains 
also substantial discussion on Z2 gauge theories and the 
doubled semion model enriched by Z2 symmetries. The 
number of phases they have obtained in cases restricted 
to local on-site symmetry action is exactly twice as ours. 
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The extra phases there can be obtained by stacking each 
of our phase, namely {TOO, TIO, Til, 500, 510, 501, 511} 
with a non-trivial Z2 SPT phase, leading altogether to 6 
distinct T phases and 8 5 phases. It would be of inter- 
est to understand possible extra phases also in the other 
constructions we have in the current paper by stacking 
them with SPT phases. 
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Appendix A: Phases in SEP{Z4 x Z4,Z2) with = 1 

In this appendix we explain how one may arrive at Eq. 
(75). Given the much experience one may gained with 
looking for condensable bosons in the other examples in 
this paper, here we shall be as brief as we can. Taking 
the scalar product of the independent bosons in the sets 
Ai,A2,A3,Bp, and Bq in Eq. (70) with the in 
Eq. (74), we can find all the relevant terms involving 
the parameters in d(p^ , which are the terms that are not 
immediately equal to (mod 27r), and are tabulated as 
follows. 



Ai 


bT2 + ch 


A2 


atl + ch 


As 


atl + dTi 


Bp 


(^ - cp)ti + bT2 + ctj. + bpTi - c'pti + c'ts 


B, 


(^ - dq)ti + bT2 + bqt3 + dTi 
+ (^ - d'q)ti + d'Ti 



(AI) 



where T2 — t2 ~ ^2^3 and = t^ — n2pi . Note that the 
coefficients in the last two rows of in the above equation 
must meet the constraints in Eqs. (70d) and (70e). Now 
that the only relevant parameters are ti, T2, t^, and T4, 
we can use a string [tiT2t3T4] to labels all possible phases 
in this case; hence there are 16 of them. Since the sets 
Ai to A3 are apparently much simpler than the infinite 
families of sets, we first see which among the 16 phases 
can have their edge modes fully gapped by condensing the 
bosons in these simpler sets. It is immediate from the first 
three rows in Eq. (Al), as long as any one of the pairs 
{T2,t3), (<i,t3), and {ti,T4) is (0,0), one can condense 
the bosons in the set whose relevant terms are turned 
off by the corresponding vanishing pair of parameters. 
Thus, the following phases have trivial edge modes as 



completely gapped. 



[tlT2t3T4] 


Condensable sets 


[0000] 


Any one 


[0001] 


Al or A2 


[0010] 


A3 


[0100] 


A2 or A3 


[0101] 


A2 


[0110] 


A3 


[1000] 


Al 


[1001] 


Al 



Now we have 8 phases left. We have two infinite fam- 
ilies of sets at our disposal. Consider the family Bp 
first, if we take p & Z + ^, i.e., half integers, the con- 
straints c'p € Z and bp ^ Ia in Eq. (70d) demands that 
c',6, c G 2Z; hence, we can let b = 2j, c' = 2k and 
assumed p = 1/2 for simplicity without losing any gener- 
ality, which renders another constraint {bpn2/2~cp) e Z 
in Eq. (70d) as j = 2m + c with m G Z. With these in 
mind, the two relevant terms of Bp become 

(^ - c)ti + 2jT2 + ct3 + 2mTi + cT^ - kti + 2kt3 

= -^ti + c{t3 + T4) (mod 2), 

where an overall tt factor is dropped. Then clearly, if 
ti = and t^ + Ti = 2, the equation above automatically 
hold, indicating that the edge modes in phases [0011] and 
[0111] can be completely gapped out by condensing any 
set Bpgz+i/2- 

Now let us turn to the family Bq. Since, as argued 
before in Section VII B, n2 = is equivalent to stacking 
two copies of toric code model that is studied in our first 
example, we focus on 712 = 1 from now on for simplicity. 
First consider q = (Ak — l)/2. A: G Z, according to Eq. 
(70e), this readily constrains that b = 2j € 27L and that 
d! 12 — d! q — 2kd! + d! . Hence, the two relevant terms of 
Bg in Eq. (Al) become 

,^ 4fc-l, 4fc-l 

(f - rf^— )*i + 2jl^2 + 2j^— 13 + 

+ {2kd! -\- d')ti + d'Ti 
= {d + d'){ti + Ti) - jt3 (mod 2) 

where 0.2 = 1 is assumed. Thus, ilt^ — and ti+T^ = 2, 
the two relevant terms will become irrelevant, implying 
that the edge modes in the phase [1101] can be fully 
gapped by condensing any Bqg2Z-i/2- 

One can verify by similar procedures that condensing 
any set Bq with q G 2Z + 1/2, the phase [1100] has a 
completely gapped edge without breaking the symmetry. 
We catalog the above results in the following table. 



[hT2t3Ti] 


Condensable sets 


[0011] 


Bp, pe2Z 


-1/2 


[0111] 


Bp, p G 2Z 


-1/2 


[1100] 


B„ g G 2Z 


f 1/2 


[1101] 


Bg, qe2Z 


f 1/2 
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We remark that the second column in above is not meant 
to be complete, in the sense that other choices of p and/or 
q may also do the job. But the point is that no set 
of independent bosons can condense without breaking 
symmetry to gap the edge modes of the remaining four 
phases: 

[tiT2hT4 = [1010], [1011], [1110], [nil]. (A4) 

Appendix B: Some useful matrices 

K matrices of the form 



for / € {0, 1, • • ■ , iV — 1} feature frequently in our discus- 
sion of topological phases which descend from deconfined 
Q = Zn gauge theories. 

We give a list of SL{2,'E) matrices X{N,l) that keep 
K{N,l) invariant . 

There arc three special cases where there are general 
solutions of X. 



[N,0] — O-x, 



1 

-N -1 



N -1 N -2 



X[N,1] = 

^[N,N-1] - \ 1 _ TV 

More generally, it can be parametrized as 

h/k {h^-k^)s/N 



(B2) 



X 



[N,l] 



-N/{k^s) -h/k 



(B3) 



where I — h x k x s for h,k, s G Z. The parameteriza- 
tion follows from Euler's paramctrization of Pythagorean 
triples. Wc note that not all I therefore allow for an 
X[N,i]eSLi2,Z). 



Appendix C: From SPT to Topological phases 

As noted first in Ref"'^° and elaborated in Ref^'*'^^, 
there is a close relation between a bosonic SPT phase 
with symmetry group Gs and a topological phases char- 
acterized by a deconfined gauge group Q where Q = Gs- 
The precise relation one can turn an SPT phase with 
only short range entanglement into a topological gauge 
theory with long range entanglement by introducing an 
extra set of gauge fields and gauging Gs- 

This procedure has a direct analog also in the context 
K matrix construction. 

Recall that a generic bosonic SPT phase can be con- 
structed by taking K = 3S> the starting point and 
then imposing global symmetry by incorporating suitable 



Higgs terms that respect the symmetry'^. The Chcrn- 
Simons Lagrangian is thus 



(CI) 



Let us be specific and consider in particular SPT phases 
with Zat symmetry. In that case, the symmetry trans- 
formation is characterized by Ref'^ 



{W^dct>^} = {t,2^/N{l,q)}. 



(C2) 



This dictates how the anyonic excitations character- 
ized by charge vector 1 created by bi ^ cxp(iZ^07) trans- 
form. Recall also that the Chern-Simons construction is 
the "dual frame" description of these bosons'^^, where the 
currents of these bosons are related to the CS gauge fields 
by ill ~ ie/^i/p/(27r)9yap. (This is a standard normaliza- 
tion. See e.g. Ref'^'^'^''.) Therefore we can write down 
the current of the Zjy symmetry in terms of a^, which is 
given by 



27r 



(C3) 



Following the standard procedure, we gauge the Zat 
symmetry by minimally coupling it to a gauge field 



^gaugo 



(C4) 



where we understand that while Aj^ is a [/ (1) gauge field, 
we are preserving only a Zn subgroup by restricting Aj^ 
to take discrete values 2TTa/N, for some integer < a < 
N- 

In the topological gauge theory, we have conservation 
of both the electric and magnetic charges. We should in- 
troduce therefore another gauge field A^ that couples to 
magnetic excitations of the "global turned local" symme- 
try. As is already evident in the discussion in Ref'''*'^^, 
the excitation of the gauge fields of the gauged Zjv is 
responsible for generating these magnetic configurations. 
Another way to see that is that in the igaugc term, by 
an integration by parts dA^ becomes electric sources of 
the a^'^ gauge fields, and it is well known that the elec- 
tric charges of a correspond to vortex excitations of the 
bosons b alluded to above in the "direct" frame. There- 
fore we introduce the following coupling 



^magnetic 



N 
An 



(C5) 



The normalization is also dictated by the fact that we 
expect unit electric charge coupled to A^ should gain a 
phase of 2t:/N when it moves around a unit magnetic 
charge coupled to A^. (c.f. discussion in Ref*). 
The total Lagrangian is then given by 



L = L 



K 



L 



magnetic- 



(C6) 



Finally, let us integrate out a^'^. Since this is a quadratic 
action, this procedure can be most readily done by ob- 
taining their equations of motion from the total action 
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and then evaluating L on-shell. The equations of motion 
are 



da^ + idA^ =0 da^ + iqdA^ = 0. (C7) 



We end up with 



L = -'-flLA^lCud^Al, (C8) 



where 

which is the expected K matrix of the topological phase 
corresponding to a (deconfined) topological gauge theory 
with gauged Zjv- 

This procedure can be readily checked for other 
bosonic SPT phases. One could readily check that the 
same procedure works for more general Abelian symme- 
try groups, such as Z2 x Z2. We note that our gauging 
procedure depends on the fact that a Z^r symmetry can 
be understood as a subgroup of U(l), which admits a nat- 
ural gauging procedure. For more general non-Abelian 
discrete symmetries we believe an analogous procedure 
should exist by embedding it in a non- Abelian Lie group. 
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